Spontaneous Pseudospin Spiral Order in the Bilayer Quantum Hall Systems 
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Using exact diagonalization of bilayer quantum Hall systems at total filling factor vt — 1 in 
the torus geometry, we show that there is a new long-range interlayer phase coherence due to 
spontaneous pseudospin spiral order at interlayer distances larger than the critical value at which the 
pseudospin ferromagnetic order is destroyed. We emphasize the distinction between the interlayer 
phase coherence and the pseudospin ferromagnetic order. 

PACS numbers: 73.43.-f, 73.21.-b 



I. INTRODUCTION 

Bilayer quantum Hall systems at total filling factor 
Vt = 1 exhibit one of the most interesting many-body 
correlation effects: spontaneous interlayer phase coher- 
ence, which is solely caused by the Coulomb interaction 
in the limit of zero interlayer tunneling. Recently this 
phenomenon has received drastic experimental support 
by Spielman et alX, who discovered a strong enhance- 
ment in the zero-bias interlayer tunneling conductance 
for interlayer distance d < d c ~ 1.83/p i n the regime of 
very little tunneling, where Iq = yJhc/eB is the magnetic 
length. These experiments strongly indicated a presence 
of spontaneous interlayer coherence at small interlayer 
distances. 

The ground state in the limit of zero interlayer distance 
can be shown to be spontaneously interlayer phase coher- 
ent due to layer symmetry. One of the most convenient 
way to see this is to use the pseudospin representation in 
which the top (bottom) layer is denoted as pseudospin 
| (|). In the limit of zero interlayer separation, there 
is exact pseudospin SU(2) symmetry so that our bilayer 
system at vt = 1 with pseudospin can be exactly mapped 
onto a single layer system v — 1 with real spin. Then, it 
is straightforward to see that the ground state should be 
a pseudospin ferromagnet in order to reduce the Coulomb 
energy cost, i.e. Hund's rule. On the other hand, it has 
been known that, in the limit of large interlayer separa- 
tion, the ground state is composed of two split composite 
fermion Fermi seas^*^. 

Therefore, it is clear that there should be a quantum 
phase transition at a critical interlayer distance d c since 
the ground states at small and large interlayer separa- 
tions are not adiabatically connected. However, it is not 
at all clear what happens near d c . In spite of intense the- 
oretical efforts in the past, the true nature of phase tran- 
sition near d c has been elusive. Analytical approaches 
based on Hatree-Fock theories&i were not able to reli- 
ably treat strong quantum fluctuations near the phase 
transition, which led to the erroneous prediction of bro- 
ken translational symmetry. On the other hand, ana- 
lytic approaches utilizing effective field theories^ were 
not directly based on the microscopic Hamiltonian of 
the many-body Coulomb interaction. So its validity 



should be justified by other means such as exact diago- 
nalization. In particular, mapping of the original many- 
body Coulomb Hamiltonian to a spin model Hamiltonian 
is not expected to be valid near the critical interlayer 
separation. Previous numerical studies based on exact 
diagonalizatioriiiSiii, however, did not take into account 
fundamental fluctuations of interlayer number difference 
due to intrinsic, quantum-mechanical uncertainty in layer 
indices. Consequently those studies could not provide 
any direct information related to spontaneous interlayer 
phase coherence. 

It is crucial at this point to distinguish between pseu- 
dospin ferromagnetic order and interlayer phase coher- 
ence. There is no fundamental reason to believe that 
they are identical because the interlayer phase coherence 
is a necessary condition, but not a sufficient condition for 
pseudospin ferromagnetism. We will show in this article 
that, at intermediate interlayer distances, there is an- 
other sort of long-range interlayer phase coherence due 
to pseudospin spiral ordering. It should be emphasized 
that this pseudospin spiral order is spontaneous in that 
it exists solely because of the Coulomb interaction with- 
out any external parallel magnetic field. The formation 
of pseudospin spiral order has been considered previously 
in the presence of parallel magnetic field£ii2. 

This article is organized as follows. We will begin by 
defining the Hamiltonian of bilayer quantum Hall sys- 
tems in Sec. ^ The pseudospin magnetization is com- 
puted as a function of interlayer distance d in the limit 
of zero interlayer tunneling in Sec. IIIII where the intrin- 
sic, quantum-mechanical fluctuation of interlayer number 
difference is fully taken into account, based on the anal- 
ogy with superconductivity. Defined as the point where 
the pseudospin ferromagnetic order is destroyed, the crit- 
ical interlayer distance d c is reliably determined. In Sec. 
IIVI we focus on the nature of the instability which de- 
stroys pseudospin ferromagnetic order at d c . We will 
show that the instability is due to the formation of pseu- 
dospin spiral order which is manifested in the lowest en- 
ergy excitations. By computing a new order parameter 
for pseudospin spiral, we will show in Sec. El that the 
ground state itself has the pseudospin spiral order. Fi- 
nally, it is argued in Section IVII that a signature of the 
phase transition from pseudospin ferromagnetic order to 
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our predicted pseudospin spiral order may have been al- 
ready observed experimentally. 



II. HAMILTONIAN 

Let us begin by writing the Hamiltonian for the bilayer 
quantum Hall systems: 



H = H t + PlllVcouiPlll 



(1) 



where Plll is the lowest Landau level (LLL) projection 
operator. Vcoul represents the usual Coulomb interaction 
between electrons: 
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(2) 

where = y/Hc/eB is the magnetic length, d is the inter- 
layer distance, and is the lateral separation between 
the i-th and j-th electrons. In the above we have used a 
pseudospin representation to distinguish the top (|) and 
the bottom (|) layers. In general, we define the pseu- 
dospin operator as follows: 



S = 



(3) 



where a is the usual Pauli matrix, and m denotes the LLL 
orbital index. Note that S z is half the electron number 
difference between the two layers, and S x is associated 
with interlayer tunneling. We take the real spin to be 
fully polarized. 

The tunneling Hamiltonian H t in Eq. Q can be written 
as: 
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where t is the single particle interlayer tunneling gap. 
Although Eq.QJ is valid for general t, we are inter- 
ested only in the limit of zero interlayer tunneling, i.e. 
t/(e 2 /elo) — > 0, which is appropriate when considering 
spontaneous interlayer coherence (note that the t — > 
limit is not the same as the t = situation). We analyze 
the Hamiltonian in Eq. Q by using exact diagonalization 
(via a modified Lanczos method) in the torus geometry^. 



III. SPONTANEOUS PSEUDOSPIN 
MAGNETIZATION 

As mentioned in the beginning, it will be shown that 
there is a new long-range interlayer phase coherence due 
to pseudospin spiral order at sufficiently large d. But, 
first, we should determine the critical interlayer distance, 
d c , at which pseudospin ferromagnetic order terminates. 
Of course, the most natural order parameter for the pseu- 
dospin ferromagnetic order is the pseudospin magnetiza- 
tion in the cc-direction: (S x ). 



It is tempting to apply exact diagonalization tech- 
niques to finite systems with fixed number of electrons 
in each layer. But, S x changes the interlayer number dif- 
ference by unity (AS Z = ±1) so that the ground state ex- 
pectation value is precisely zero when computed naively 
without any explicit interlayer tunneling. This problem 
has been addressed in an ad hoc manner by introducing 
explicit interlayer tunneling 11 , which, however, severely 
obscures the effect of spontaneous phase coherence and 
may produce misleading results. The real solution is to 
realize that there is intrinsic, quantum-mechanical uncer- 
tainty in the layer index of electrons at small interlayer 
distance (even in the limit of zero tunneling) so that the 
true ground state is a linear combination of various 
states with different S z : 
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where \<Pm) is the lowest energy state with S z = M, Am 
is a sharply peaked function of M with width 0(y/~N), 
and N is the number of electrons. It is important to re- 
member that even the Halperin (1,1,1) state, ^(1,1,1), has 
Am oc exp (— 2M 2 /N), which can be proved by analyzing 
the following, explicit wavefunction including both the 
orbital and pseudospin part 15 : 



i^(i,i,i)> =n 
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^ (1,1,1) is the exact ground state at d — 0. Because it 
can be exactly mapped onto the BCS wavefunction, it is 
fruitful to consider an analogy with the BCS theory for 
general d. Following the BCS theory, one can compute 
the pseudospin magnetization S x as follows: 



(lp\S x \lp) = ^ ^M^M>((f>M>\S x \(f>M} 



M,M' 

— (4>M*+l\S x \(t>M*) + ((f>M*-l\S x \(t>M*) 

- 2(4>M*+l\S x \4>M*) (7) 

where the last two step are justified in the thermody- 
namic limit since Am is sharply peaked at M — M*. 
Although (S x ) should be exactly N/2 at d = in the 
thermodynamic limit due to the exact pseudospin SU(2) 
symmety, there are some finite-system size corrections. 
After some algebra on the (1,1,1) state, Eq.Q shows that 
(^(1,1,1)1 ^x1^(1,1,1)) is (iV + l)/2forNodd (M* = —1/2), 
and is y/N(N + 2)/2 for N even (A/* = 0). Therefore, 
the pseudospin magnetization should be scaled as N + 1 
(^iV(7V + 2)) for N odd (even). 

FigH shows the pseudospin magnetization per parti- 
cle (defined in Eq.Q) as a function of d/l n for various 
numbers of electrons. As conventional in finite system 
analysis, the critical interlayer distance can be obtained 
from the inflection point of (S x ): d c ~ 1.5?o- Although 
d c was estimated to be around 1.5/o previously in various 
approaches, it should be emphasized that the pseudospin 
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magnetization (S x ) is computed for the first time in this 
article without any interlayer tunneling. And therefore 
d c is reliably determined without any ambiguity. 

IV. PSEUDOSPIN SPIRAL INSTABILITY 

We now address the nature of instability causing the 
destruction of pseudospin ferromagnetic order. This 
question is, in turn, directly related to the nature of the 
new ground state at d > d c . To be specific, we con- 
sider the lowest energy excitations which are responsible 
for the ground state instability. Fig|2| shows the disper- 
sion curves of the lowest energy excitations at d/lo = 0.5 
and 1.5 for a finite system with the total number of elec- 
trons N = 13 and the number of pseudospin-up electrons 
= 7 (or —6). There are several points to be em- 
phasized. 

First, there is clearly a linearly dispersing Goldstone 
mode at small interlayer distances, such as d/lo = 
0.5. However, withinnumericalaccuracyvegaxdviig the 
discreteness of finite-system wavevectors, the Goldstone 
mode seems to vanish for d > d c , which coincides with the 
destruction of pseudospin ferromagnetic order—. Consis- 
tent with the conclusion from (S x ), this shows that the 
ground state undergoes a phase transition around d = d c . 

Second, contrary to the prediction of the random phase 
approximation (RPA) theory 6 , there is no softening of 
the dispersion curve (roton) at any finite momentum for 
any interlayer distance. Remember that the RPA theory 
predicts a roton around kl$ ~ 1.4, which was believed 
to cause an instability toward charge density wave order. 
This already shows that previous theories have some se- 
rious flaws in describing the phase transition at d = d c . 
We will show in the next paragraph that there are com- 
pletely different low-energy excitations at large d/lo. 

Third, several previous numerical studies^iii were un- 
fortunately based on the exact diagonalization of N — 
8 system, which, under careful investigation^, can be 
shown to predict a Wigner crystal at large d rather than 



FIG. 2: Dispersion curves of the lowest energy excitations 
at interlayer distance d/lo = 0.5 and 1.5 for a finite system 
with the total number of electrons N — 13 and the number 
of pseudospin-up electrons AT-j = 7 (or Afy =6). It will be 
shown that the lowest energy state in the momentum channel 
k/o = y2ft/N(Ni , 0) (denoted by arrows) is the pseudospin 
spiral state. 



a quantum Hall liquid, which is a finite-system artifact 
of boundary condition and special number of electrons. 

Fourth, these lowest-excitation energies are exactly 
zero in the d — > oo limit where there is no interaction 
between electrons with different pseudospins. 

However, the most crucial aspect of the dispersion 
curve is that there arc completely new low-energy ex- 
citations due to the pseudospin spiral order, which are 
denoted by arrows in Fig|5] The lowest energy state in 
the momentum channel kl = y / 2n/N (JVj, , 0) (denoted 
by arrows in Fig|2J) begins to be ripped out of the well- 
defined dispersion curve at d/lo — 0.5, and it becomes 
completely separated from all the other excitations at 
d/lo > 1.5 so that it becomes really the lowest energy 
excitation among all momentum channels. The same be- 
havior is found in all studied systems with different values 
of N and , with exception of the N = 8 system due to 
a finite-size artifact as mentioned in the above. Note that 
there are three other degenerate excitations in the mo- 
mentum channels kZ = v / 2Tr/N(N h 0), ■ s /2ir/N(Q,N i ), 
and ^2tt/N(Q, iVj) due to reflection symmetries and pe- 
riodic boundary conditions of torus geometry. These pe- 
culiar excitations are completely different from usual ro- 
tons because (i) they form isolated points rather than 
a part of conventional dispersion curve, and (ii) their 
momentum increases as |k|Zo cx \/N, which diverges in 
the thermodynamic limit. We will show later that these 
peculiar properties are precisely the properties of pseu- 
dospin spiral state. 

In order to understand why these excitations are the 
pseudospin spiral states, it would be best if we first ex- 
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FIG. 3: Overlap between the pseudospin spiral wavefunction 
(described in the text) and the lowest energy state in the 
Hilbert space of momentum channel kZo = ^2tt/N(Ni, 0). 



amine the following pseudospin spiral wavefunction: 



\ipspiral(n)) = Yl 



V2 ( " 



"ml 



10) (8) 



where m denotes the linear momentum in the Landau 
gauge, i.e., p x = 2irm/L and L is the linear system size 
(note that L = 1^2-kN at vt = 1). We will show 
later in FigOU by explicitly computing the overlap be- 
tween \ip sp i ra i(n = 1)) and the lowest-energy excitations 
(denoted by arrows in FigEJ) that the lowest-energy ex- 
citations are the pseudospin spiral states. But, first, we 
would like to emphasize that the pseudospin spiral state 
I'fpspiraiin)) is the exact ground state in the presence of 
parallel magnetic field 9 . A surprising result of this ar- 
ticle is that, even without parallel magnetic field, spiral 
states are the lowest energy excitations which cause the 
instability of pseudospin ferromagnetic order. 

To elucidate the physical meaning of pseudospin spiral 
state, examine Eq.JSJ which reveals that \ip spiral) con- 
tains diagonal interlayer correlations between the states 
with p x — 2n(m + n)/L and 2i:m/L. Remember that 
the interlayer correlation between p x = 2i:{m + n)/L 
and 2-Km/ L is identical to the diagonal interlayer correla- 
tion between electrons with different pseudospins which 
are separated laterally in the y-direction by rd§\J 2n j N . 
In turn, this diagonal interlayer correlation is physically 
equivalent to the pseudospin spiral order with spiraling 
period of L/n because the pseudospin part of wavefunc- 
tion is essentially given by 
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(9) 



Pseudospin spiral excitations (and therefore the diago- 
nal interlayer correlation) are schematically depicted in 
FigH 

Now we compute the momentum of pseudospin spi- 
ral states by directly applying the magnetic translation 
operator onto \ip spiral ("■))• The explicit algebra shows 
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FIG. 4: Schematic diagram for (a) the pesudospin ferromag- 
netic state (note the vertical interlayer correlation), (b) the 
pseudospin spiral state with clockwise spiraling direction, and 
(c) the pseudospin spiral state with anti-clockwise spiraling 
direction. Ellipses indicate the Coulomb correlation of neu- 
tral pairs between electron in one layer and hole in the other 
layer. It is shown in the text that there is an interlayer co- 
herence due to pseudospin spiral order at d > d c , where the 
ground state is conjectured to be a bound state of two pseu- 
dospin spiral states with opposite winding direction, i.e. a 
bound state of (b) and (c). 



that the momentum of ip S pirai{n) is |k|Zo = niVf 2tt /N 
which is exactly the momentum of the lowest-energy ex- 
citation seen in FigJ^Jif n = 1. Note that the momentum 
of this pseudospin spiral state is actually proportional 
to its spiraling period. That is the reason why the mo- 
mentum of pseudospin spiral state with long spiraling 
period diverges in the thermodynamic limit. Similarly, 
we can explain why the pseudospin spiral state occurs 
in an isolated momentum channel. Consider two pseu- 
dospin spiral states with long period, say n = 1 and 
n = 2, whose momenta differ by a factor proportional to 
\/~N /Iq. Therefore, the pseudospin spiral states with two 
different long spiraling periods do not occur in adjacent 
momenta. 

It is important to note that \tl> s piral{n = 1)) is the pseu- 
dospin spiral excitation in the limit of long spiraling pe- 
riod (= L), which initiates the destruction of pseudospin 
ferromagnetic order. Therefore, it may be indicative of a 
second order phase transitioni 9 -. 

As mentioned earlier, Fig0 shows the overlap be- 
tween the pseudospin spiral state \ip S pi r ai{n = 1)) in 
Eq.© and the lowest-energy excited state with k/o = 
^2tt/N(N^, 0)2!i. As one can see, the overlap is high 
overall at small d/lo, and is very close to unity (~ 0.97) 
especially around d/lo — 0.5 which coincides with the 
point where (S x ) begins to deviate from the fully sat- 
urated value of 1/2. This shows that the instability of 
pseudospin ferromagnetic order is indeed caused by pseu- 
dospin spiral order. 
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FIG. 5: Ground state expectation value of the new order 
parameter (T„ with n = 1) measuring the diagonal interlayer 
correlation. Note that the diagonal interlayer correlation is 
physically equivalent to the pseudospin spiral order. 

It should be noted that, though the overlap is lowered 
at larger d/lo, the low energy excitation is in general adi- 
abatically connected to the pseudospin spiral state. In 
fact, the reason for lower overlaps at larger d/lo is that 
our trial wavefunction in Eq. ||SJ is not very accurate any- 
more for larger d/lo since it is based on the assumption 
that the ground state is the (1,1,1) state. Of course, 
the ground state itself is affected by the pseudospin spi- 
ral instability which, after all, destroys the pseudospin 
ferromagnetic order at the critical interlayer distance. 
Therefore, it would be satisfactory if, in addition to the 
evidence due to the low energy excitation, one can prove 
directly that the ground state has the pseudospin spiral 
order at d ~ d c . The next section will be devoted to the 
pseudospin spiral order in the ground state. 

V. PSEUDOSPIN SPIRAL ORDER IN THE 
GROUND STATE 

It is conjectured that the ground state at large d is in 
fact a bound state of two pseudospin spiral states with 
opposite winding direction, i.e. a bound state of (b) and 
(c) in Fig^J which requires a complicated interaction be- 
tween the two spiral states (note that this conjecture is 
consistent with the fact that the ground state itself does 
not break translational symmetry). So instead of con- 
structing a trial wavefunction we take a more transpar- 
ent approach by defining a new order parameter which 
measures the pseudospin spiral order (or equivalently the 
diagonal interlayer correlation): 

Tn = 2 51 ( C L+n,T C ^i + 4ni c m+n,]) (10) 
m 

The following properties of T n are worth mention- 
ing, (i) T n is a generalization of pseudospin magneti- 
zation S x which measures the vertical interlayer correla- 



tion: T n —o = S x - So the relationship between T n and 
ipspirai{n) is analogous to that between S x and t/'fi.i.i) • 

(ii) When acting upon momentum eigenstates, T n al- 
ters momentum of the eigenstates, though changes in mo- 
mentum are negligible in the thermodynamic limit. So, 
in order to avoid any confusion, we compute (T 2 ) instead. 

(iii) In order to be meaningful in the thermodynamic 
limit, (T 2 ) should be scaled as iV 2 . 

(iv) Since the (1,1,1) state does not have any 
pseudospin spiral order, ("0(1,1,1) l^nlV'(i,i,i))/-W' 2 
should be zero in the thermodynamic limit. How- 
ever, there are some finite-system size corrections: 

|22Mi,i,i)> = N/4 - N 1 N l /2(N - 1) = O(N). 
We subtract these finite system size corrections from the 
expectation value of T 2 . 

FigEl plots the ground state expectation value of T 2 =1 
as a function of interlayer distance d/lo- As one can see, 
(T% =1 )/N 2 is peaked around d/lo — 1.5 which is con- 
sistent with the critical interlayer distance d c . So it is 
shown that indeed the ground state contains the pseu- 
dospin spiral order at d ~ d c . In fact, FigEl implies that 
the pseudospin ferromagnetic order and the pseudospin 
spiral order may coexist for a range of d since the ex- 
pectation value of pseudospin spiral order parameter is 
not zero even in interlayer distances smaller than d c . Of 
course, in the thermodynamic limit, (T 2 =1 )/iV 2 will be 
zero at d — ► oo since there is no interlayer correlation in 
the two split Fermi seas of composite fermions. 

Finally, we would like to discuss the physical origin of 
pseudospin spiral order. First, let us remind ourselves 
that the pseudospin ferromagnetic state can be viewed 
as a condensate of neutral pairs of electrons and corre- 
lation holes which are bound directly across the inter- 
layer separation. Now for sufficiently large d the ground 
state has the pseudospin spiral order, or equivalently di- 
agonal interlayer correlation, which can be also viewed 
as a condensate of neutral electron-hole pairs with ex- 
tended size and higher symmetry such as p-wave, as com- 
pared to the s-wave-like pair in pseudospin ferromagnet. 
Now, the origin of spiral order can be identified with the 
origin of extended pairs which is due to weak pairing. 
Of course, the pairing between electrons and correlation 
holes is weakened as d increases. Eventually the bind- 
ing is completely lost when the split composite fermion 
Fermi seas are formed. 



VI. CONCLUSIONS 

We conclude by mentioning some experimental impli- 
cations. Recently, Kellogg et al. observed an enhanced 
longitudinal drag resistance at d/lo &s large as 2.6 where 
the enhanced tunneling conductance is absent 21 . Theo- 
retically, tunneling is renormalized from the bare single 
particle tunneling gap t to the greatly enhanced renor- 
malized tunneling t(S x ). Therefore, it is natural to as- 
sume that the enhanced tunneling conductance becomes 
absent as soon as the pseudospin ferromagnetic order is 
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destroyed at d = d c . On the other hand, we have shown 
that there is a new interlayer phase coherence due to the 
pseudospin spiral order for d > d c which, we speculate, 
causes the longitudinal drag anomaly because of remain- 
ing interlayer correlation. Regarding the Hall resistance, 
our pseudospin spiral state is expected to be incompress- 
ible, as shown in the dispersion curve. However, the low- 
est energy gap is estimated to be very small so that the 



Hall plateau will be difficult to be observed in current 
temperature ranges and impurity concentrations. 
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